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Abstract. We review recent results on global attractors of U(l)-invariant dispersive Hamil- 
' Ionian systems. We study several models based on the Klein-Gordon equation and sketch 

, the proof that in these models, under certain generic assumptions, the weak global attractor 

is represented by the set of all solitary waves. In general, the attractors may also contain 
multifrequency solitary waves; we give examples of systems which contain such solutions. 
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Q ■ 1 Introduction 

, 

The long time asymptotics for nonlinear wave equations have been the subject of intensive 
^ ■ research, starting with the pioneering papers by Segal [Ml [45], Strauss [53], and Morawetz and 

. Strauss [39], where the nonlinear scattering and local attraction to zero were considered. Global 

attraction (for large initial data) to zero may not hold if there are quasistationary solitary wave 

1 J ■ solutions of the form 

■ ij{x,t) = ^x)e''''^\ with u; G M, lim (/>(a;) = 0. (1.1) 

We will call such solutions solitary waves. Other appropriate names are nonlinear eigenf unctions 
and quantum stationary states (the solution (II. ip is not exactly stationary, but certain observable 
quantities, such as the charge and current densities, are time-independent indeed). 

Existence of such solitary waves was addressed by Strauss in [53] , and then the orbital stability 
of solitary waves has been studied in [211 SS SHI US] • The asymptotic stability of solitary waves 
has been studied by Sof fer and Weinstein [501 IS] i Buslaev and Perelman [3 [8] , and then by 
others. 

The existing results suggest that the set of orbitally stable solitary waves typically forms 
a local attractor, that is to say, attracts any finite energy solutions that were initially close to it. 
Moreover, a natural hypothesis is that the set for all solitary waves forms a global attractor of 
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all finite energy solutions. We address this question in the present paper, reviewing the results 
on the global attraction in several models based on the Klein-Gordon equation, and describing 
the developed techniques. 

We briefly discuss the long-time solitary wave asymptotics for U(l)-invariant Hamiltonian 
systems in Section [2j The definitions and results on global attraction to solitary waves from the 
recent papers [MlEllET] are presented in Section [3l We also give there a very brief sketch of the 
proof. In Section HI we give a description of all the steps (omitting excessive technical points) 
of the argument for the simplest model: Klein-Gordon equation interacting with a nonlinear 
oscillator. The key parts of the proof are presented in full detail. The examples of (untypical) 
multifrequency solitary waves are given in Section [5j 



2 History of solitary asymptotics 

for U(l)-invariant Hamiltonian systems 

2.1 Quantum theory 

Bohr's stationary orbits as sohtary vi^aves 

Let us focus on the behavior of the electron in the Hydrogen atom. According to Bohr's postu- 
lates [U, an unperturbed electron runs forever along certain stationary orbit, which we denote \E) 
and call quantum stationary state. Once in such a state, the electron has a fixed value of ener- 
gy E, not losing the energy via emitting radiation. The electron can jump from one quantum 
stationary state to another, 

\E^) ^ \E+), (2.1) 

emitting or absorbing a quantum of light with the energy equal to the difference of the ener- 
gies E^ and E^. The old quantum theory was based on the quantization condition 

p • dq = l-Khn, n G N. (2.2) 

This condition leads to the values 

4 

me 

for the energy levels in Hydrogen, in a good agreement with the experiment. Apparently, the 
condition (|2.2|) did not explain the perpetual circular motion of the electron. According to 
the classical Electrodynamics, such a motion would be accompanied by the loss of energy via 
radiation. 

In terms of the wavelength A = of de Broglie's phase waves [6], the condition ()2.2p states 
that the length of the classical orbit of the electron is the integer multiple of A. Following de 
Broglie's ideas, Schrodinger identified Bohr's stationary orbits, or quantum stationary states \E), 
with the wave functions that have the form 

V^(x,t) =,/.^(x)e-*"*, io = E/h, 

where h is Planck's constant. Physically, the charge and current densities 

— e — — 

p{yi,t) = eipip, j(x,t) = — (^/j • - VV' • V'), 

with e < being the charge of the electron, which correspond to the (quasi)stationary states of 
the form ip{'x.,t) = (/)tj(x)e~*^* do not depend on time, and therefore the generated electromag- 
netic field is also stationary and does not carry the energy away from the system, allowing the 
electron cloud to flow forever around the nucleus. 
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Bohr's transitions as global attraction to solitary waves 

Bohr's second postulate states that the electrons can jump from one quantum stationary state 
(Bohr's stationary orbit) to another. This postulate suggests the dynamical interpretation of 
Bohr's transitions as long-time attraction 

^{t) — > \E±), t ±00 (2.3) 

for any trajectory ^{t) of the corresponding dynamical system, where the limiting states \E±) 
generally depend on the trajectory. Then the quantum stationary states Sq should be viewed as 
the points of the global attractor £/. 

The attraction (j2.3p takes the form of the long-time asymptotics 

il^{x, t) ~ (x)e-^'^±*, t ^ ±00, (2.4) 

that hold for each finite energy solution. However, because of the superposition principle, the 
asymptotics of type (j2.4p are generally impossible for the linear autonomous equation, be it the 
Schrodinger equation 

^2 e^ 
ihdtip = -—Alp - —ip, 
2m |x| 

or relativistic Schrodinger or Dirac equation in the Coulomb field. An adequate description 
of this process requires to consider the equation for the electron wave function (Schrodinger 
or Dirac equation) coupled to the Maxwell system which governs the time evolution of the 
four-potential A{x,t) = {ip{x,t), A{x,t)): 

{ihdt — e(p)'^fp = ( c-V — eA j ip + rn^c^ilj, 

Dip = A7re{tpi;-5{x)), DA = 47re ^ " " " ^ . 

Consideration of such a system seems inevitable, because, again by Bohr's postulates, the tran- 
sitions (|2.ip are followed by electromagnetic radiation responsible for the atomic spectra. More- 
over, the Lamb shift (the energy of 2S1/2 state being slightly higher than the energy of 2P^/2 
state) can not be explained in terms of the linear Dirac equation in the external Coulomb field. 
Its theoretical explanation within the Quantum Electrodynamics takes into account the higher 
order interactions of the electron wave function with the electromagnetic field, referred to as 
the vacuum polarization and the electron self-energy correction. 

The coupled Maxwell-Schrodinger system was initially introduced in [43]. It is a U(l)- 
invariant nonlinear Hamiltonian system. Its global well-posedness was considered in [22]. One 
might expect the following generalization of asymptotics (12. 4p for solutions to the coupled 
Maxwell-Schrodinger (or Maxwell-Dirac) equations: 

{^{x, t),A{x, t)) ~ {^^^ (x)e-'"±*, A^^ (x)) , t^±oo. (2.5) 

The asymptotics (|2.5p would mean that the set of all solitary waves 

{{(P^e-'^^^A^) -.LoeR} 

forms a global attractor for the coupled system. The asymptotics of this form are not available 
yet in the context of coupled systems. Let us mention that the existence of the solitary waves 
for the coupled Maxwell-Dirac equations was established in [T6] . 
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2.2 Solitary waves as global attractors for dispersive systems 

Convergence to a global attractor is well known for dissipative systems, like Navier-Stokes equa- 
tions (see [U [531 [SB]), such systems, the global attractor is formed by the static stationary 
states, and the corresponding asymptotics (|2.4p only hold for t — > +00. 

We would like to know whether dispersive Hamiltonian systems could, in the same spirit, 
possess finite dimensional global attractors, and whether such attractors are formed by the 
solitary waves. Although there is no dissipation per se, we expect that the attraction is caused 
by certain friction mechanism via the dispersion (local energy decay) . Because of the difficulties 
posed by the system of interacting Maxwell and Dirac (or Schrodinger) fields (and, in particular, 
absence of the a priori estimates for such systems), we will work with simpler models that share 
certain key properties of the coupled Maxwell-Dirac or Maxwell-Schrodinger systems. Let us 
try to single out these key features: 

(z) The system is \J{1) -invariant. 

This invariance leads to the existence of solitary wave solutions 4>^{x)e~'^^^ . 

(a) The linear part of the system has a dispersive character. 

This property provides certain dissipative features in a Hamiltonian system, due to local 
energy decay via the dispersion mechanism. 

(Hi) The system is nonlinear. 

The nonlinearity is needed for the convergence to a single state of the form 4'^^{x)e~'^^^ . 
Bohr type transitions to pure eigenstates of the energy operator are impossible in a linear 
system because of the superposition principle. 

We suggest that these are the very features are responsible for the global attraction, such as (|2.4p . 
()2.5p , to "quantum stationary states" . 

Remark 1. The global attraction (j2.4p . (12. 5p for U(l)-invariant equations suggests the corre- 
sponding extension to general G-invariant equations (G being the Lie group): 

^(x, t) ~ ^±(x, t) = e"±V±(2;), t ±00, (2.6) 

where Q.± belong to the corresponding Lie algebra and e^±* are the one-parameter subgroups. 
Respectively, the global attractor would consist of the solitary waves (j2.6p . On a seemingly 
related note, let us mention that according to Gell-Mann-Ne'eman theory [18] there is a cor- 
respondence between the Lie algebras and the classification of the elementary particles which 
are the "quantum stationary states". The correspondence has been confirmed experimentally 
by the discovery of the omega-minus Hyperon. 

Besides Maxwell-Dirac system, naturally, there are various nonlinear systems under consid- 
eration in the Quantum Physics. One of the simpler nonlinear models is the nonlinear Klein- 
Gordon equation which takes its origin from the articles by Schiff [41^ I42| . in his research 
on the classical nonlinear meson theory of nuclear forces. The mathematical analysis of this 
equation is started by Jorgens [26] and Segal [1^, who studied its global well-posedness in the 
energy space. Since then, this equation (alongside with the nonlinear Schrodinger equation) has 
been the main playground for developing tools to handle more general nonlinear Hamiltonian 
systems. The nonlinear Klein-Gordon equation is a natural candidate for exhibiting solitary 
asymptotics ()2.4p . 

Now let us describe the existing results on attractors in the context of dispersive Hamiltonian 
systems. 
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Local and global attraction to zero 

The asymptotics of type (j2.4p were discovered first with 'ip± = in the scattering theory. Namely, 
Segal, Morawetz, and Strauss studied the (nonlinear) scattering for solutions of nonlinear Klein- 
Gordon equation in |39l B6l [53] . We may interpret these results as local (referring to small 
initial data) attraction to zero: 

ijj{x, t) ~ V'i = 0, ±oo. (2.7) 

The asymptotics ()2.7p hold on an arbitrary compact set and represent the well-known local 
energy decay. These results were further extended in [19l [201 113 EH ■ Apparently, there could be 
no global attraction to zero {global referring to arbitrary initial data) if there are solitary wave 
solutions 0a;(x)e~*'^*. 

Solitary waves 

The existence of solitary wave solutions of the form 

with i/^(]R") being the Sobolev space, to the nonlinear Klein-Gordon equation (and nonlinear 
Schrodinger equation) in M"-, in a rather generic situation, was established in [53] (a more general 
result was obtained in [21 [3] ) . Typically, such solutions exist for to from an interval or a collection 
of intervals of the real line. We denote the set of all solitary waves by Sq. 

Due to the U(l)-invariance of the equations, the factor-space 5o/U(l) in a generic situation 
is isomorphic to a finite union of intervals. Let us mention that there are numerous results on 
the existence of solitary wave solutions to nonlinear Hamiltonian systems with U(l) symmetry. 
See e.g. [UlinillZ]. 

While all localized stationary solutions to the nonlinear wave equations in spatial dimensions 
n > 3 turn out to be unstable (the result known as "Derrick's theorem" [E]), quasistationary 
solitary waves can be orbitally stable. Stability of solitary waves takes its origin from [58] and 
has been extensively studied by Strauss and his school in [21[ Wl\ [^ [l9] . 

Local attraction to solitary waves 

First results on the asymptotics of type (j2.4p with uj± ^ were obtained for nonlinear U(l)- 
invariant Schrodinger equations in the context of asymptotic stability. This establishes asymp- 
totics of type (j2.4p but only for solutions close to the solitary waves, proving the existence of 
a local attractor. This was first done by Softer and Weinstein and by Buslaev and Perelman 
in [3 [HI [501 [51], and then developed in [HI [IH [H [H [IS [52] and other papers. 

Global attraction to solitary waves 

The global attraction of type ()2.4p with ifj-j- ^ and lv± = was established in certain models 
in [281 [29l [30l \3T\ [32] [33] for a number of nonlinear wave problems. There the attractor is the 
set of all static stationary states. Let us mention that this set could be infinite and contain 
continuous components. 

In [34] and [35], the attraction to the set of solitary waves (see Fig. [T|) is proved for the 
Klein-Gordon field coupled to a nonlinear oscillator. In [36], this result has been generalized 
for the Klein-Gordon field coupled to several oscillators. The paper [37] gives the extension to 
the higher-dimensional setting for a model with the nonlinear self-interaction of the mean field 
type. We are going to describe these results in this survey. 
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Figure 1. For t ±00, a finite energy solution '^{t) approaches the global attractor £/ which coincides 
with the set of all solitary waves Sq. 

We are aware of but one recent advance [55] in the field of nontrivial (nonzero) global attrac- 
tors for Hamiltonian PDEs. In that paper, the global attraction for the nonlinear Schrodinger 
equation in dimensions n > 5 was considered. The dispersive (outgoing) wave was explicitly 
specified using the rapid decay of local energy in higher dimensions. The global attractor was 
proved to be compact, but it was neither identified with the set of solitary waves nor was proved 
to be finite-dimensional \55\ Remark 1.18]. 

3 Assumptions and results 

In [351 ESI [37] we introduce the models which possess the key properties we mentioned above: 
U(l)-invariance, dispersive character, and the nonlinearity. The models allow to prove the global 
attraction to solitary waves and to develop certain techniques which we hope will allow us to 
approach more general systems. 

Model 1: Klein— Gordon field with a nonlinear oscillator 

We consider the Cauchy problem for the Klein-Gordon equation with the nonlinearity concen- 
trated at the origin: 

ii}{x, t) = Tp"{x, t) - m^^{x, t) + (5(x)F(^/>(0, t)), x G M, 

ip\t=o = ^0(2;), ^L^o = (3.1) 

Above, ?n > and F is a function describing an oscillator at the point x = 0. The dots stand 
for the derivatives in t, and the primes for the derivatives in x. All derivatives and the equation 
are understood in the sense of distributions. We assume that equation (j3.ip is U(l)-invariant; 
that is, 

F(e^V) = e^V('0), 6* G M, V G C. (3.2) 

If we identify a complex number ip = u + iv ^ C with the two-dimensional vector (n, v) G M^, 
then, physically, equation (j3.ip describes small crosswise oscillations of the infinite string in 
three-dimensional space (x, u, v) stretched along the x-axis. The string is subject to the action 
of an "elastic force" —m?"ilj{x,t) and coupled to a nonlinear oscillator of force Fiyijj) attached at 
the point x = 0. We assume that the oscillator force F admits a real-valued potential, 

F(V^) = -V^f/(^), G C, [/ G C2(C), 

where the gradient is taken with respect to Re'i/' and Imil^. 

Remark 2. Viewing the model as an infinite string in M^, the assumption (13. 2p means that the 
potential U (ip) is rotation-invariant with respect to the x-axis. 
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The codes for the numerical simulation of a finite portion of a string coupled to an oscillator, 
with transparent boundary conditions {dfip = idxi^), ^-re available in 

Model 2: Klein— Gordon field with several nonlinear oscillators 

More generally, we consider the Cauchy problem for the Klein-Gordon equation with the non- 
linearity concentrated at the points Xi < X2 < ■ ■ ■ < X^: 

N 

i>{x, t) = t) - m^i;{x, t)+Y, - Xj)Fj{i;{Xj,t)), x G M, 

J=i 

V-lt^o = ^0(2;), tpl^o = Mx)- (3.3) 

Model 3: Klein— Gordon field with the mean field interaction 

We also consider the Klein-Gordon equation with the mean field interaction: 

ii}{x,t) = Aip{x,t) - m^ij{x,t) + p{x)F{{p,^p{-,t))), a; G M", n > 1, 

ipl^o = M^), V'lt.o = ^0(2;). (3.4) 

Above, /9 is a smooth real- valued coupling function from the Schwartz class: p G o5^(]R"), p ^ 0, 
and 



(p,^(-,0) = / /o(x)V(x,t) 



dx. 



(3.5) 



Hamiltonian structure 

Equations (|3.4p formally can be written as a Hamiltonian system, 

where ^ = {ip,Tr) and DTC is the Frechet derivative of the Hamilton functionals 

Wosc(V',vr) = i j {\7r\^ + \i^f+m^\i;\^)dx + Y.^jiijiXj)), 

R 

Wn,.f.(V',vr) = ^ I {\7r\^ + \Vi;\^ + m^\i;\')dx + U{{p,i;)). 

Since (j3.3p and (13. 4p are U(l)-invariant, the Nother theorem formally implies that the value of 
the charge functional 

Q{'4>^ ~ 2 j ~ 

is conserved for solutions ^{t) = (■0(i), vr(t)) to (j3.5p . 

Let us introduce the phase space (f of finite energy states for equations (|3.3p , (|3.4p . Denote 
by II ■ ||l2 the norm in the complex Hilbert space L^(IR") and by || • ||^2 the norm in L^(B^), 
where is a ball of radius R > 0. Denote by || • ||hj the norm in the Sobolev space //"(B^) 
(which is the dual to the Sobolev space //(^**(B^) of functions supported in the ball of radius R). 
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Definition 1 (The phase space). 

(i) S = H'^{W) © l2(]R"), n > 1, is the Hilbert space of the states with the norm 

||(V^,^)|||:= \\m, + \\V^\\l,+m^i:\\l,. 

(ii) For e > 0, (f"" = H^-^{W) ® iJ-'(M") is the space with the norm 
,n)U-. = \\{l-Ar^/\^|;,n)U. 



[iii] Define the seminorms 



?p ^ is the space with the norm 

oo 

. = 5^2-^||(V^,vr)|U-.,^. (3.6) 



oo 

^)ll<r-^ = 

R=\ 



Equations (|3.3|) . (|3.4p are formally Hamiltonian systems with the Hamilton functionals "Hose 
and ?Ym.f.) respectively, and with the phase space S from Definition [1] (for equation (j3.3p . the 
dimension is n = 1). Both Hose (or Wm.f.) and Q are continuous functionals on S . 

Global well-posedness 



Theorem 1 (Global well-posedness). Assume that the nonlinearity in (|3.3p is given by Fj{z) 
—Wj{z) with inf L 

respectively). Then: 



VUj{z) with inf Uj{z) > -oo, 1 < J < N (or F{z) = -VU{z) with inf U{z) > -oo in 



(z) For every (■iPqjTTq) G S' the Cauchy problem (|3.3p ( (|3.4p . respectively) has a unique global 
solution 'il){t) such that (^,^) G C{K,S). 

(ii) The map W{t) : {■iPqjTTq) i— > {ip{t),-ip{t)) is continuous in S for each t G IR. 

(iii) The energy and charge are conserved: 7i{ip{t),'ip{t)) = const, Q{ijj{t),ip(t)) = const, t G M. 

(ii;) The following a priori bound holds: 

\\{m,m)h<C{i'o,TTo), tGM. (3.7) 
(v) For any e G [0, 1], 

(iA,iA)gc(^)(m,#-^), 

where C^'''^ denotes the space of Holder- continuous functions. 



The proof is contained in [351 ES] , and [37] . 
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Solitary waves 

Definition 2 (Solitary waves). 

(z) The solitary waves of equation (j3.3p are solutions of the form 

= 0^(x)e-*'^*, where w G M, 0^ G F^(M"). (3.8) 

(ii) The set of ah solitary waves is So = {(puj- uj ^M., (puj ^ H^{W^)]. 
(Hi) The solitary manifold is the set S = {{(t)un —i^'Puj)' u; G M, (pui ^ //^(M")} C S". 

Remarlc 3. 

(i) The U(l) invariance of (j3.3p and (|3.4p implies that the sets Sq, S are invariant under 
multiplication by e*^, 9 £ M. 

(ii) Let us note that for any G M there is a zero solitary wave with (j)uj{x) = since Fj{0) = 0. 

The following proposition provides a concise description of all solitary wave solutions to (j3.ip . 

Proposition 1. There are no nonzero solitary waves for \lo\ >m. 

For a particular uj G {—m,m), there is a nonzero solitary wave solution to (\'S.l\] if and only 
if there exists C G C\{0} so that 

2k{uj) = F{C)/C. (3.9) 
The solitary wave solution is given by 

= C7e-''('^)l^l, K{io) = -u;2. (3.10) 
Remarlc 4. There could be more than one value C > satisfying ()3.9p . 
Remark 5. By (|3.10p . lo = itm can not correspond to a nonzero solitary wave. 
Proof. When we substitute the ansatz (^c^e"*'^* into (13. Ih . we get the following relation: 

-uj^Mx)=(l)'u(.x)-m^^^{x)+5{x)F{Mx)), x G M. (3.11) 

The phase factor e"*"^* has been canceled out. Equation p. lip implies that away from the origin 
we have 

hence (f>uj{x) = C±e~'^^^^^ for ±x > 0, where k± satisfy = — uj"^. Since (f>uj{x) G H^, it is 
imperative that k± > 0; we conclude that \uj\ < m and that k± = \J vn? — uP- > 0. Moreover, 
since the function (j)io{x) is continuous, C_ = C+ = C ^ (since we are looking for nonzero 
solitary waves). We see that 

(^^(x) = Ce-^l^l, K = - lo^ > 0. (3.12) 

Equation (13. lip implies the following jump condition at a; = 0: 

= ci)Uo+) - <PUo-) + F{MO)), 

which is satisfied due to p.9p and ()3.12p . ■ 
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Global attraction to solitary waves 

We will combine the results for the Klein-Gordon equation with one and several oscillators ( (|3.1 
and (1331)). 



Theorem 2 (Global attraction for (jS.Sp . Klein-Gordon equation with N oscillators). Assume 
that all the oscillators are strictly nonlinear: for all 1 < J < N , 

FjOj) = -VUjiij), (3.13) 
PJ 

where C/j(V') = ^ uj,i^^, uj,pj > 0, and pj>2. 

1=0 

Further, if N > 2, assume that the intervals [Xj, Xj^i], 1 < J < N — 1, are small enough so 
that 

/ ^2 \ 1/2 /J N \ 

where pj are exponentials from (j3.13|) . Then for any {i^QjITq) £ the solution ^(t) to the 
Cauchy problem (13. 3p converges to S: 



lim dist^^((V(t),^(t)),5) =0. (3.15) 

t— ►±oo 

Theorem 3 (Global attraction for (j3.4p . Klein-Gordon equation with mean field interaction). 
Assume that the nonlinearity F{z) is strictly nonlinear: 

p 

F{z) = -VU{z), where U{z) = ^ui\z\'^\ u/ e M, Up > 0, and p>2. 

1=0 

Further, assume that the set 

Zp = {u; £ R\[-m, m]: p(C) = for ah ^ G M" such that m'^ + f = lv^} 
is finite and that 

does not vanish at the points uo G Zp. Then for any (^O)^o) £ <^ solution 'ip{t) G C{M,£') to 
the Cauchy problem (j3.4p converges to S in the space (^f'', for any e > 0; 

^lim dist^^.(4'(t),5) = 0. (3.16) 

Above, dist «-£ (^', 5) := inf 11^' — $|L-e, with 11 • |L-£ is defined in (13.611 . 

Theorem [2] is proved in [36]; Theorem [3] is proved in [37]. We present the sketch of the proof 
of Theorem [2] for one oscillator in Section [H 
Let us mention several important points. 

(i) In the linear case, the global attractor contains the linear span of points of the solitary 
manifold, (S). In [35], we prove that for the model of one linear oscillator attached to the 
Klein-Gordon field the global attractor indeed coincides with (S). 
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(ii) The condition (j3.14p allows to avoid "trapped modes", which could also be characterized 
as multifrequency solitary waves. In Proposition [S] below, we give an example of such 
solutions in the situation when the condition (j3.14p is violated. 

Similarly, the condition of Theorem [3] that cr{uj) does not vanish for uj & Zp allows to avoid 
multifrequency solitary waves. 

(iii) We prove the attraction of any finite energy solution to the solitary manifold S: 

{m,m)^s, t^±oc, (3.17) 

where the convergence holds in local seminorms. In this sense, 5 is a weak (convergence 
is local in space) global (convergence holds for arbitrary initial data) attractor. 

(iv) S can be at most a weak attractor because we need to keep forgetting about the outgoing 
dispersive waves, so that the dispersion plays the role of friction. A strong attractor would 
have to consist of the direct sum of S and the space of outgoing waves. 

(v) We interpret the local energy decay caused by dispersion as a certain friction effect in 
order to clarify the cause of the convergence to the attractor in a Hamiltonian model. 
This "friction" does not contradict the time reversibility: if the system develops backwards 
in time, one observes the same local energy decay which leads to the convergence to the 
attractor as t — > — oo. 

(vi) Although we proved the attraction ()3.17|) to S, we have not proved the attraction to 
a particular solitary wave, falling short of proving (12.41) . Hypothetically, if 5/U(l) contains 
continuous components, a solution can be drifting along S, keeping asymptotically close to 
it, but never approaching a particular solitary wave. This could be viewed as the adiabatic 
modulation of solitary wave parameters. Apparently, if 5/U(l) is discrete, a solution 
converges to a particular solitary wave. 

(vii) The requirement that the nonlinearity is polynomial allows us to apply the Titchmarsh 
convolution theorem. This step is vital in our approach. We do not know whether the 
polynomiality requirement could be dropped. 

(viii) For the real initial data, we obtain a real-valued solution 'ip{t). Therefore, the conver- 
gence (IXT5D . (I3T6]1 of ^{t) = {'i/j{t),ip{t)) to the set of pairs (0^, -iujcj)^) with lo G M\{0} 
implies that 'ip{t) locally converges to zero. 



Sketch of the proof 

First, we introduce a concept of the omega- limit trajectory f3{x,t) which plays a crucial role in 
the proof. 



Definition 3 (Omega-limit trajectory). The function P{x,t) is an omega-limit trajectory if 
there 

that 



there is a global solution ^ G C(M, £') and a sequence of times {sj: j G N} with lim Sj = oo so 



^(x,t -I- Sj) i3{x,t), 
where the convergence is in Ch{[—T, T] x B^) for any T > and R > 0. 

We are going to prove that all omega-limit trajectories are solitary waves: /3(x, t)=0(^(x)e~*'^*. 
It suffices to prove that the time spectrum of any omega-limit trajectory /? consists of at most 
one frequency. 

To complete this program, we study the time spectrum of solutions, that is, their complex 
Fourier-Laplace transform in time. First, we prove that the spectral density of a solution is 
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absolutely continuous for > m hence the corresponding component of the solution disperses 
completely. It follows that the time-spectrum of omega-limit trajectory (3 is contained in a finite 
interval [—m,m]. 

Second, we notice that /? also satisfies the original nonlinear equation. Since the spectral 
support of (3 is compact and the nonlinearity is polynomial, we may apply the Titchmarsh con- 
volution theorem. This theorem allows to conclude that the spectral support of the nonlinearity 
would be strictly larger than the spectral support of the linear terms in the equation (which 
would be a contradiction!) except in the case when the spectrum of the omega-limit trajectory 
consists of a single frequency G [—m,m]. 

Since any omega-limit trajectory is a solitary wave, the attraction (j3.17p follows. 

Open problems 

(z) As we mentioned, we prove the attraction to S, as stated in (|3.17p . but have not proved the 
attraction to a particular solitary wave like (j2.4p . It would be interesting to find solutions 
with multiple omega-limit points, that is, the situation when the frequency parameter lv 
keeps changing adiabatically. 

(ii) Our argument does not apply to the Schrodinger equation. The important feature of the 
Klein-Gordon equation is that the continuous spectrum corresponds to |a;| > m, hence the 
spectral density of the solution is absolutely continuous for \uj\ > m, while the spectrum 
of the omega-limit trajectory is within the compact set [— m, m]. This is not so for the 
Schrodinger equation: since the continuous spectrum corresponds to w > 0, the resulting 
restriction on the spectrum of the omega-limit trajectory is u; < 0. As a result, we do 
not know whether the spectrum is compact; the Titchmarsh convolution theorem does not 
apply, and the proof breaks down. It would be extremely interesting to investigate whether 
the convergence to solitary waves is no longer true, or instead certain modification of the 
Titchmarsh theorem allows to reduce the spectrum to a point. 

(iii) Similarly, the Titchmarsh theorem does not apply when the nonlinearity is not polynomial, 
and it would be interesting to investigate what could happen in such a case. 

4 Proof of attraction to solitary waves 

for the Klein— Gordon field with one nonlinear oscillator 

We will sketch the proof of Theorem [2] for the system (|3.ip which describes one nonlinear 
oscillator located at the origin. 

Proposition 2 (Compactness. Existence of omega-limit trajectories). 

(i) For any sequence sj — > -|-cxd there exists a subsequence Sji +oo such that 



for some /3 G C(M x M), where the convergence is in Cb{[—T, T] x [—R, R]), for any T > 



tp{x, Sji 



+ t)^ /3{x,t) 



(4.1) 



and R > 0. 






Proof. By Theorem [T] (w), for any e G [0, 1], 
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Taking e = 1/4, we see that ip £ C^"^(]R x M), for any a < 1/4. Now the first statement of 
the proposition follows by the Ascoli-Arzela theorem. The bound (j4.2|) follows from (|4.1|) . the 
bound (j3.7p . and the Fatou lemma. ■ 

We call omega-limit trajectory any function (3{x, t) that can appear as a limit in (j4.ip (cf. 
Definition [3]) . We are going to prove that every omega-limit trajectory /? belongs to the set of 
solitary waves; that is, 

= (/)^_^(x)e~*'^+* for some uj+e[-m,m]. (4.3) 

Remark 6. The fact that any omega-limit trajectory turns out to be a solitary wave implies 
the following statement: 

if there is a sequence tj ^ oo so that e X L'^, then $ G 5. (4.4) 

In turn, (j4.4p implies the convergence to the attractor in the metric (j3.6p of S'p'' for e > 0: 

{7p{t),iP{t)) ^S, e > 0, t^ ±oo. 

This is weaker than the convergence to the attractor in the topology of £'p stated in Theorem [21 
For the proof of the convergence to the attractor in the topology of (op, see |35j . 



Let us split the solution 'il){x,t) into two components, il){x,t) = xix,t) + '^{.x,t), which are 
defined for all t G M as solutions to the following Cauchy problems: 

X{x,t) = x"{x,t) - m^x(.x,t), {x,x)\t=o = (^o(2;),vro(x)), 

^{x,t) = ip"{x,t)-m^ip{x,t) + 6{x)f{t), ((^,(^)|^^„ = (0,0), (4.5) 

where {ip(){x),TrQ{x)) is the initial data from ()3.ip . and 

fit) := F{i;{0,t)), teR. (4.6) 

Lemma 1 (Local energy decay of the dispersive component). There is the following decay for 
X- 

lim ||(x(-,0,x(-,i))IUi? = 0, Vi?>0. 

For the proof, see [351 Lemma 3.1]. Lemma [1] means that the dispersive component x does 
not give any contribution to the omega- limit trajectories (see Definition [3]) . 



^(cj — zO) = + -\/ cj^ —mp- ■ ^ , k{uj — iO) = — ^/uj^ — m-^ 

fc(ti;itiO)=ivm?— <^ 



Figure 2. Domain D and the values of k{uj ± iO), a; e M. 
Let k{uj) be the analytic function with the domain D := C\((— oo, —m] U [m, +oo)) such that 
k{uj) = \I<J^ - m'^, lmfc(u;)>0, to £ D. (4.7) 
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See Fig. [2j Let us also denote the limit of k{u}) for u + iO, a; G M, by 

k+{oj) := k{oo + iO), w G M. (4.8) 
The function f{x, t) = ip{x, t) — xi^7 1) satisfies the following Cauchy problem: 

(p{x, t) = ip"{x, t) - rr?ip{x, t) + 6{x)f{t), {if, ip)\^^,, = (0, 0), 

with f{t) defined in (j4.6p . Note that ■0(0, •) G Cb(M) by the Sobolev embedding, since (ipjip) G 
Cb(M, (#) by Theorem [1] (ill). Hence, f{t) G Cb(M). On the other hand, since xi^^'^) is a finite 
energy solution to the free Klein-Gordon equation, (x, x) ^ Ch(^, S')- It follows that = — x 
is also of finite energy norm: 

(<^,vi) G C6(M,^). (4.9) 

We denote 

(^+(x,t) := e{t)^{x,t), U{t) := e{t)f{t) = 0(t)(V'(O,t)). 
The function Lp^{x,t) satisfies the equation 

y5+{x,t) = dlip+{x,t)-m^^+{x,t) + 6{x)U{t), {^+,ip+)l^, = (0,0), teR. (4.10) 

We set Tt^^[g{t)]{uj) = / e''^^g{t)dt for a function g{t) fr om the Schwartz space o$^(M). The 
Fourier transform 



oo 



^+{x,uj) = J^t->ujW+{x,t)] = / e''^'^{x,t)dt, {x,oj) G M^ 

JO 

is a continuous function of x G M with values in tempered distributions of w G M, which satisfies 
the following equation (cf. (j4.10p ): 

—(jo'^(fj^{x, uj) = d'^(fj^{x, Lo) — ?n^(^+(x, Lo) + 6{x)f+{io), (x, oj) G M?. 
Proposition 3 (Spectral representation). There is the following relation: 

(^+(x,a;) = (^+(0,w)e*'=+('^)l^l, x G M. (4.11) 
Proof. Let us analyze the complex Fourier transform of (/j+(x,t): 

POO 

<f+{x,uj) = J^t~.uj[ip+{x,t)]= e''^^ip{x,t)dt, xGM, G C+, 



where C"*" := {z G C : Imz > 0}. Due to (j4.9p . (^+(-,w) are ff^-valued analytic functions of 
UJ G C"*". Equation (j4.10p implies that (f+ satisfies 

— a;^(^+(x, Lv) = 9^(^+(x, io) — m^Lp^{x, iv) + (5(x)/+(a;), u> G C"*". 

^±ik{uj)\x\ 

The fundamental solutions G±(x,uj) = — — r- satisfy 

±2ik[uJ) 

Gl{x,uj) + {uj^ -m^)G±{x,uj) = 6{x), uj G C+. 

Note that for each oj G C"*" the function G+(-,u;) is in H^{M.) by definition (j4.7p . while G^{-^uj) 
is not. The solution (^_|_(x,a;) can be written as a linear combination of these fundamental 
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solutions. We use the standard "limiting absorption principle" for the selection of the appropriate 
fundamental solution: Since (p+{-,uj) € H^(K) for u £ C"*", so is G+(-,a;), while G_(-,a;) is not, 
we have: 

ip+{x,u;) = -U{u;)G+{x,u;) = ^, u; G C+. (4.12) 

The relation KW\ yields 

(^+(x,u;) = -/+H^-— - = e^'=('^)l^l<^+(0,u;), x G M, u; G C+. (4.13) 

Now we extend the relation (Ii33|) to a; G M. Since if G Cb(R,H^(R)) by (gj]), we have 
e{t)^{x,t)= lim e{t)ip{x,t)e-^\ (4.14) 

where the convergence holds in the space of i7^-valued tempered distributions, =5^'(R, iif^(M)). 
The Fourier transform ip^{x,uj) := Tt-tui[(p+{x,t)] = J^t~*oj[d{'t)^{x,t)] is defined as a tempered 
If^-valued distribution of a; G M. As follows from (j4.14p and the continuity of the Fourier 
transform J-'t^co in =5^'(M), ip+{x,uj) is the boundary value of the analytic function (^+(x,u;), in 
the following sense: 

(f+{x,uj)= lim ip+{x,uj + ie) = lira J^t^^[9{t)(p{x,t)e~'^^], w G M. (4.15) 

Again, the convergence is in the space =5^'(M, i/^(M)). 

We use (j4.15p to take the limit Ima; — > 0+ in the expression (j4.13|) for (^+(x,u;), and keep 
in mind that (p+{x,uj) is a quasimeasure (see Remark [7]) for each x G M, while the exponential 
factor in (j4.13p is a multiplicator in the space of quasimeasures. The formula (j4.1ip follows. ■ 

Remark 7. A tempered distribution /x(a;) G is called a quasimeasure if fi{t)=,^~\^^[fi{uj)] 

G Cb(M). For more details on quasimeasures and multiplicators in the space of quasimeasures, 
see [351 Appendix B]. 

Proposition 4 (Absolute continuity of the spectrum). The distribution (p^{0,uj) is absolutely 
continuous for \uj\ > m, and moreover 



j 



(^+(0,u;)|^ ' duj < oo, (4.16) 



where k+{uj)/uj > for to e M\[-?n, m] (see (|4.8p and Fig.^. 

Proof. We use the Paley-Wiener arguments. Namely, the Parseval identity and (j4.9p imply 
that 

oo 

j \\^+i-,u; + ie)\\l2duj = 27r j e-^'^\\ip+{-,t)\\l2dt<j, e > 0. (4.17) 

R 

On the other hand, we can calculate the term in the left-hand side of ()4.17p exactly. According 

to Km . 

<^+{x,uj + ie) = (^+(0,u; + ie)e^'=('^+^")l^l, 
hence (|4.17p results in 

e / |(^+(0,cu + ie)|2||e^'=(^+*")l^l||2,dc^<C7, e > 0. (4.18) 
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Here is a crucial observation about the norm of e*'^('^+*^)l^l. 
Lemma 2. 

(i) For u) E 



lim e||e^'=(^+^^)l^l||f2 

e->0+ 



/c+(a;) 



(m) For any 5 > there exists £,5 > such that for \u>\ > m + 6 and e G (0, e^), 

k+{u:) 



^||gifc(a;+je)|a;| ||2^ > 



2UJ 



(4.19) 



(4.20) 



Remark 8. The asymptotic behavior of the L^-norm of e*'^('^+*^) stated in the lemma is easy to 
understand: for lo € IR\[— m, m], this norm is finite for e > due to the small positive imaginary 
part of k{LO + ie), but it becomes unboundedly large when e — > 0+. Let us also mention that 
the integral (j4.19p is easy to evaluate in the momentum space. 



Substituting (OH into (fiJ8]) . we get: 
/ |(^+(0,cj + ie)|2.^±Mdc^ < 2C, 

J llx) 



UJ 



< e < es, 



(4.21) 



\ui\>m+5 

with the same C as in (j4.18p . We conclude that for each 6 > the set of functions 



gs,si^^) = (^+(0,0; + ie) 



UJ 



1/2 



defined for u G i?^, is bounded in the Hilbert space L^(]R\[— m — 5,m + 5]), and, by the Banach 
Theorem, is weakly compact. The convergence of the distributions (j4.15p implies the following 
weak convergence in the Hilbert space L^(R\[— m — 5,m + S]): 



95,1 



95, 



0+, 



where the limit function gs{Lo) coincides with the distribution (^+(0, w) 
M\[— m — 5,m + 6]. It remains to note that, by (I4.2ip . the norms of al 



1/2 



restricted onto 



functions gs, 5 > 0, are 



bounded in 



-m — 6,m + 6]) by a constant independent on 6, hence ()4.16p follows. 



By Lemma [H the dispersive component xi'^t) converges to zero in S'p as t ^ 00. On the 
other hand, by (j4.ip . 7p{x,t + Sj') converges to P{x,t) as j' — > 00, uniformly on every compact 
set of the plane M^. Hence, ip{x, t + Sj') = Tp{x, t + Sj/) — x{x, t + Sji) also converges to (3{x, t), 
uniformly in every compact set of the plane M^: 



(p{x,Sj> +t) ^ l3{x,t), 



x G 



t G 



(4.22) 



Therefore, taking the limit in equation (14. Sp . we conclude that the omega-limit trajectory /3(x, t) 
also satisfies the same equation: 



(3{x, t) = (3"{x, t) - m^f3{x, t) + 6{x)F{(3), x G 



t G 



Taking the Fourier transform of /? in time, we see by (14. ip that [3{x,lo) is a continuous function 
of X G M, with values in tempered distributions of a; G M, and that it satisfies the corresponding 
stationary equation 



-00'^ ${x,uj) = P"{x,uj) — m'^P{x,uj) + 6{x)g{u!), {x,uj) G 



(4.23) 
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valid in the sense of tempered distributions of {x, lo) G M^, where g{uj) are the Fourier transforms 
of the function 

git) :=F(/3(0,t)). 
For brevity, we denote 

m :=/5(0,t). 
Lemma 3 (Boundedness of spectrum). 

supp (3 C [—m, m]. 
Proof. By (g^S), we have 

ip+{x, Sj> +t) ^ I3{x, t), X £R, t£ M, (4.24) 
with the same convergence as in (j4.ip and (j4.22p . We have: 

(f+ {x,Sj +t) = — I e-*'^*e"*'^^^' (p+{x,u})dui, x G M, t G M, 

where the integral is understood as the pairing of a smooth function (oscillating exponent) with 
a compactly supported distribution. Hence, (14.24p implies that 

e"*'^V(^+(x,u;) ^ /3(x,u;), x£R, Sj, ^ oo, (4.25) 

in the sense of quasimeasures (the convergence in the space of quasimeasures is equivalent to 
the Ascoli-Arzela type convergence of corresponding Fourier transforms; see |35t Appendix B]). 
Since (^+(0, cj) is locally for \uj\ > m by Proposition [H the convergence ()4.25p at x = shows 
that (3(u;) := /3(0,u;) vanishes for \uj\ > m. This proves the lemma. ■ 

We denote 

k{u) := -ik+{uj), wGM, (4.26) 
where k-^.{uj) was introduced in ()4.8p . We then have ReK^uj) > 0, and also 

k{u;) = \/ io"^ — m? > for — m < uj < m, 
in accordance with (j3.10p . 

Proposition 5 (Spectral representation for /5). The distribution P{x,uj) admits the following 
representation: 

/3(x,u;) = |3(w)e-'^('^)l^l, x G R. 

Proof. This follows by taking the limit in the first line of (j4.1ip . since supp (3 C [—m,m] by 
LemmaO while k{u>) = iHi{uj) for —m < u! < m (cf. (j4.26p ). ■ 

Proposition 6 (Reduction to point spectrum). Either supp |3 = for some uj-^- G [—m,m] 

or |3 = 0. 
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Proof. By Lemma [3l we know that supp (3 C [—m,m]. According to equation (|4.23p . the 
function (3 satisfies the foUowing jump condition at the point x = 0: 



Since supp/3'(0ib, •) C supp (3 by Proposition O it follows that 



suppff(-) C supp p. 



(4.27) 



On the other hand, by (|3.13|) . the Fourier transform g^uj) of g{t) := F{P{0,t)) is given by 



Now we will use the Titchmarsh convolution theorem [57] (see also |38^ p. 119] and [24:\ 
Theorem 4.3.3]) which could be stated as follows: 

For any compactly supported distributions u and v, supsupp(ti * v) = sup supp n + sup supp z;. 

Applying the Titchmarsh convolution theorem to the convolutions in (j4.28p . we obtain the 
following equality: 

sup supp 5 > sup supp |3 + (p — 1) (sup supp (3 — inf supp P), (4.29) 

where we used the relation sup supp |3 = — inf supp (3. We wrote ">" because of possible cancel- 
lations in the summation in the right-hand side of ()4.28p . Note that the Titchmarsh theorem is 
applicable to (|4.28p since supp (3 is compact by Lemma [3l 
Comparing (|4.27p with ()4.29p . we conclude that 

{p — 1) (sup supp |3 — inf supp |3) = 0. 

Since p > 2 by (j3.13p (which means that the oscillator at x = is nonlinear), we conclude that 
supp |3 consists of at most a single point C [— m,m]. ■ 

By Proposition [6l supp |3 C with a;+ 6 [— m, m]. Therefore, 

|3(w) = ai(5(u; — 0;+), with some oi G C. (4.30) 

Note that the derivatives (5^*^^ (a; — oJ-y.), /c > 1 do not enter the expression for |3(tij) since (3(t) = 
/3(0, t) is a bounded continuous function of t due to the bound ()4.2p . Proposition [5] and (|4.30|) 
imply that the omega-limit trajectory /3(x, t) is a solitary wave: 

I3{x,t) = ^{x)e-'^+\ 

where (/> e F^(M) by ([42]) . This completes the proof of (jilS]) . 

Remark 9. ujj^ = itm could only correspond to the zero solution by Remark [5j 



p 




n-l 



(4.28) 



Global Attraction to Solitary Waves 



19 



5 Multifrequency solitons 

5.1 Linear degeneration 

Let us consider equation ()3.3p with N = 2, under condition (j3.14p . 

Proposition 7. If in (j3.13p one has pj = 1 for some J, then the conclusion of Theorem^ may 
no longer be correct. 

Proof. We are going to construct the multifrequency solitary waves. Consider the equation 

^ = ^" - m^^p + (5(x)Fi(^) + d{x - L)F2{'4j), 
where 

Fi(V') =ail) + /^IV'I V, ^2 W = 7-0, a, A 7 e ^■ 

Note that the function F2 is linear, failing to satisfy (j3.13p (where one now has p2 = 1). The 
function 



il){x,t) 



' (^ + B)e''('^)^sinu;t, x < 0, 
^j^^-k{u,)x ^^Kiuj)x^ ujt + C sinh(ft;(3a;)x) sin 3ujt, xe[0,L], 

C 

/^g-«(a;) _^ q^k{uj){2L-x)\ g^^^^^ _^ ^g-/.(3a;)(x-L) gin3a;i, X>L, 

' sinh(K(3a;)L) 



where uj G (0, m/?,), will be a solution if the jump conditions are satisfied at x = and at x = L: 

-ip'{0+, t) + ^'(0-, t) = aV(0, t) + /3V^(0, t), (5.1) 
-^'{L+, t) + il)'{L-,t) = aiPiL, t) + /?'03(L, t). (5.2) 

Using the identity 

o 3 1 

sin^6l = -sin6' - -sin36l, (5.3) 
4 4 

we see that 

a{A + B) siniot + /3((A + B) sinujtf 

a(A + B) + 13^^^^-^) sin cut - sin3a;t. 



4 

Collecting the terms at sinwt and at sinScot, we write the condition (|5.ip as the following system 
of equations: 

2K(a;)^ = fa(^ + S) + /3^^^^^^V (5.4) 



(A -I- R"*^ 

-k(3l<j)C = -P ^ ^ ^ . (5.5) 



Similarly, the condition ()5.2p is equivalent to the following two equations: 

2SAt(w)e'^(^)^ = 7(^e-'^('^)^ + Be"''-'^^^), (5.6) 

«;(3w)(:; 



sinh(«;(3ti;)L) 



+ K(3a;)Ccosh(K(3a;)L) = 7Csinh(K(3w)L). (5.7) 



Equations (j5.4p . (j5.5p . (j5.6p . and (j5.7p could be satisfied for arbitrary L > 0. Namely, for any 
LV G (0, m/3), one uses (15. 7p to determine 7. For any /? 7^ 0, there is always a solution A, and 
to the nonlinear system (|5.4p . (|5.6p . Finally, C is obtained from (|5.5p . ■ 
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5.2 Wide gaps 

Let us consider equation (j3.3p with N = 2. Assume that ()3.13p is satisfied. 

Proposition 8. If the condition (j3.14p is violated, then the conclusion of Theorem [2] may no 
longer he correct. 

Proof. We will show that if L := X2 — Xi is sufficiently large, then one can take Fi^ip) 
and F2{ip) satisfying (j3.13p such that the global attractor of the equation contains the multifre- 
quency solutions which do not converge to solitary waves of the form (|3.8p . For our convenience, 
we assume that Xi = 0, X2 = L. We consider the model ()3.3p with 



Fi{'tp) = F2{^^J) = F{^p), where F{tp) = a^p + ^3\^p\'^'tp, a,(3eR. 
In terms of the condition (j3.13p . pi = P2 = 2. We take L to be large enough: 

i > P.B) 
Consider the function 

= A(e-''('^)l^l +e-''('^)l^-^l)sinwt + BX[o,L](2;)sin(A:(3w)x) sinSujt, A, B eC. 
Then ilj{x^t) solves (|3.3p for x away from the points Xj. We require that 

k{2,uj) = J, (5.9) 

so that il){x,t) is continuous in rc G M and symmetric with respect io x = L/2: 

We need |u;| < m to have k{uj) > 0, and 3|ll'| > m to have k^Suj) G M. We take uj > 0, and thus 
m < 3a; < 3m. By (j5.9p . this means that we need 



7r2 



m < Y + 1^'^ < 3m. 

The second inequality is satisfied by (15. 8|] . 

Due to the symmetry of ijj{x,t) with respect to x = L/2, the jump condition both at x = 
Xi = and at x = X2 = L takes the following identical form: 

2Ak{uj) sinujt - Bk{3uj) sinSujt = F{A{1 + e""^"^^^) sinwt). (5.10) 

We use the following relation which follows from (15. 3p : 

F[A{1 + e-'^^'^)^) sinwt) = {aA{l + e-"^^^^) + ^f3\A\'^A{l + e-'^H^)3^ sinut 

- ^P\A\'^A{1 + e-^^'^^^f sin 3ujt. 
Collecting in (j5.10p the terms at sinwt and at sin3ti;t, we obtain the following system: 
2Ak{uj) = aA{l + e-^^'^^^) + ^P\A\'^A{1 + e-^(^)^)^ 

Bk{3uj) = ^l3\AfA{l + e-''('^)^)3. (5.11) 
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Assuming that ^ 7^ 0, we divide the first equation by A: 




The condition for the existence of a solution A ^ is 




(5.12) 



Once we found A, the second equation in (jS.lip can be used to express B in terms of A. 



Remark 10. Condition ()5.12p shows that we can choose (3 < taking large a > 0. The 
corresponding potential U{ip) = — a|^p/2 — /3|?/'|^/4 satisfies (|3.13p . 
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